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ABSTRACT 

We present our results of the temporal and spectral analysis of a sample of 52 
bright and hard gamma-ray bursts (GRBs) observed with the Fermi Gamma-ray 
Burst Monitor (GEM) during its first year of operation ( July 2008 - July 2009). 
Our sample was selected from a total of 253 GBM GRBs based on the event 
peak count rate measured between 0.2 and 40 MeV. The final sample comprised 
34 long and 18 short GRBs. These numbers show that the GBM sample contains 
a much larger fraction of short GRBs, than the CGRO/BATSE data set, which 
we explain as the result of our (different) selection criteria and the improved GBM 
trigger algorithms, which favor collection of short, bright GRBs over BATSE. A 
first by-product of our selection methodology is the determination of a detection 
threshold from the GBM data alone, above which GRBs most likely will be 
detected in the MeV/GeV range with the Large Area Telescope (LAT) onboard 
Fermi. This predictor will be very useful for future multiwavelength GRB follow 
ups with ground and space based observatories. Further we have estimated the 
burst durations up to 10 MeV and for the first time expanded the duration- 
energy relationship in the GRB light curves to high energies. We confirm that 
GRB durations decline with energy as a power law with index approximately 
—0.4, as was found earlier with the BATSE data and we also notice evidence of a 
possible cutoff or break at higher energies. Finally, wc performed time-integrated 
spectral analysis of all 52 bursts and compared their spectral parameters with 
those obtained with the larger data sample of the BATSE data. We find that the 
two parameter data sets are similar and confirm that short GRBs are in general 
harder than longer ones. 

Subject headings: Methods: data analysis — Gamma-ray burst: general 
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Introduction 



The most exciting results in the last decade of gamma-ray burst (GRB) science at 
MeV-GeV energies ca me from the combine d observations of the Burst And Transient Source 
Experiment (BATSEj_ 



Telescope (EGRET, 



Fishman et al. 



Fichtel et al.l 



9931 ) with the Energetic Gamma-Ray Experiment 
1994 ). Both instruments operated between 1991 and 2000 
on board the Compton Gamma-Ray Observatory (CGRO) and covered the energy bands 
from 2 0keV to 20MeV and from 20 MeV to 30 GeV, respectively. BATSE observed 2704 



bursts (jPaciesas et al 



19991 ). providing the largest GRB database from a single experiment 



thus far. Out of this sample, only fi ve bursts were detect ed with EGRET above 100 MeV 



and only one of these, GRB 930131 f lSommer et al 



1994j ). had high-energy emission that 



was consistent with an extrapolation from its spectrum obtained with BATSE between 



25 keV and 4 MeV (see also 



Dingus 



20031 ). However, later analysis of the combined 
data from the BATSE/Large Area Detectors (LADs) and the EGRET calorimeter, the 
Total Absorption Shower Counter (TASC), uncovered an MeV component in GRB 941017 



described by a 



30wer l aw with photori index approximately —1 up to about 200 MeV 



( iGonzalez et al. 



20031). 



Gonzalez et al. 



(120091 ) subsequently searched the TASC data for a 



response to 68 bright BATSE bursts. They found that only 21 showed emission detectable 
by TASC and of these, only three contained spectra with peak energy, E^ p ^v > 2 MeV. As 



Gonzalez et al. 



t00% claimed 



these spectra were only found in a time-resolved analysis, 
that the high energy component could be hidden in the brightness of the synchrotron 
emission in a time-integrated spectrum. They suggested that the existence of the high 
energy component indicated additional non-thermal processes at the source. 

With the successful launch on 2008 June 11 of the Fermi Gamma-ray Space Telescope 
{FGST, hereafter Fermi), it is now possible to search and confirm the spectral signatures 
from GRBs up to very high energies. Fermi is an international and multi-agency space 
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observatory that studies the Cosmos over an unprecedented broad energy range (lOkeV to 
300 GeV). The Gamma-Ray Burst Monitor (GBM) is the secondary instrument onboard 
Fermi, operating between ~8 keV and ~40 MeV. The search for higher energy GRB 
emission is carried out with the primary instrument, the Large Area Telescope (LAT, 



AtwoodetaL 



20091 ). a pair conversion telescope, like EGRET, operating in the energy 



range between 20 MeV and 300 GeV. GBM alone pr ovides a much wide r energy coverage 



than any other current GRB mission, such as Swift ( IGehrels et al. 



2004h . As of November 



2010, joint high-e nergy GBM/ 



100 MeV (see e.g. 



Pelassa et al 



AT o bservations have been carried out for 19 GRBs above 
20111). This number represents ~3% of the total number of 



GRBs observed with GBM. 

The small number of LAT GRB detections may be due to instrumental bias. The 
combined GBM/LAT GRB spectra are usually well-described in the MeV-Ge V range by a 



single 



power law with an index in the approximate range of —1 to —3 (see e.g. 



AbdoetaL 



2009aJ). This result is in agreement with the distributions of the high -energy power-law 



2006 



indices observed with BATSE in the ~30 keV-~2 MeV energy range (iKaneko et al. 
hereafter K06). Consequently, photon counts above ~1 MeV are usually very low, and this 
combined with the limited Field-of-View (FoV) of the LAT, results in much fewer GRBs 
observed in the multi-MeV band than in the keV-band. Currently, detailed calculations 
of LAT upper limi ts are being performed for the brightest bursts detected with GBM 



flAbdo et al. 



2011ah 



In this paper we establish a well-defined sample of 52 bright GBM GRBs with 
broad-band spectral coverage and with statistically significant high spectral resolution. All 
events were collected during the first year of GBM operations (July, 14, 2008 - July 15, 
2009) and were detected up to MeV energies with GBM alone. In Section |2] we discuss 
the GBM instrumentation and data types and in Section [3] we describe the selection 
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methodology for our sample, mainly based on the peak count rates measured above 500 
keV. In Section m we present the temporal analysis of our sample over several energy bands. 
In Section [S] we describe the results of our time-integrated spectral analysis using different 
photon models to fit each spectrum, and discuss the distributions and correlations of the 
spectral parameters for the best models. Our conclusions are presented in Section [6l 

2. Instrumentation and Data Types 

The primary role of GBM was to augment the science return from Fermi in the study of 
GRBs by making observations at lower energies (~8 keV to ~40 MeV) and thus bridging the 
gap with those of the LAT. The GBM flight hardware comprises a set of 12 Thallium-doped 
Sodium Iodide crystals (Nal(Tl), hereafter Nal) and two Bismuth Germanate crystals 
(Bi4Ge30i3, commonly abbreviated as BGO). The individual Nal detectors are mounted 
around the spacecraft in four groups of three. Their arrangement results in an exposure 
of the whole sky unocculted by the Earth in orbit. The Nal detectors are able to detect 
7-rays in the energy range between ~8 keV and ~1 MeV. The two BGO detectors are 
mounted on opposite sides of the Fermi spacecraft. With their energy range between ~0.2 
and ~40 MeV, they provide the overlap in energy with the LAT instrument and are crucial 
in the study of high-energy, hard bursts. 

To trigger the GBM flight software (FSW), two or more Nal detectors must have 
a simultaneous statistically significant rate increase above the background rate (usually 
> 5cr). This requirement increases the threshold against statistical fluctuations and 
suppresses triggering due to non-astrophysical events that appear in only one detector, such 
as phosphorescence spikes. Before performing any spectral analysis, the detector geometry 
with respect to the GRB direction must be carefully taken into account. Detectors which 
see the burst at an angle >50°, or which suffer from blockages (by the solar panels, by the 



- 6 - 



LAT or by the spacecraft itself) were excluded. Sometimes these criteria result in a single 
Nal detector to be chosen for the spectral analysis, which is then fitted together with the 
mostly illuminated BGO detector. The best detector combination for each burst is given in 
columns 4 (Nals) and 5 (BGO) of Table [H 

All three GEM data types, namely CSPEC, CTIME and TTE data, were used for the 



analysis pre sented in the 



be found in 



Meegan et al 



lowing sections. A detailed description of these data types can 



fl2009r ). CTIME and CSPEC data were used in burst-mode, 
i. e. at 64 ms and 1.024 s temporal resolution, respectively. The BGO peak count-rate 
analysis and the determination of the burst durations in the integrated BGO energy range 
are based on CTIME files, which have the finest temporal resolution (64 ms) with modest 
energy resolution consisting of 8 energy channels. For the determination of the duration 
dependence on energy, we used TTE data both for the Nal and the BGO detectors. CSPEC 
and TTE data provide an energy resolution consisting of 128 energy channels and were 
used for all spectral analyses. The coarse time-resolution CSPEC data are normally used as 
pre-trigger background data for the TTE data, since the latter only include ~ 30 s before 
the trigger time. TTE data are then collected up to 300 s post-trigger, and in all cases 
discussed hereafter cover the whole burst duration. For each GRB, the data type used for 
spectral analysis is listed in column 7 of Table [H 



3. Sample Selection Methodology 

We selected our sample based on two requirements: (i) a significant count rate excess 
above background (> 3a) measured by the most illuminated Nal detectors in the 50-300 
keV energy range, to ensure good statistics; and (ii) a significant count rate excess above 
background measured by the most illuminated BGO detector in the 500 keV-1 MeV during 
the main burst emission episode (Tgo). This combination allows a broadband spectral 
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analysis of GBM data, which spans about three decades in energy (10 keV-10 MeV). It also 
represents a significant difference from previous selection criteria adopted by K06, which 
were mainly based on peak photon flux and fluence values measured in the 50-300 keV 
energy range. Criterion (ii) was afforded by the good performance of the BGO detectors at 
those energies I > 40 keV), where the effective area of the Nal detectors rapidly decreases 
( Bissaldi et al.l 120091) . Below, we discuss the methodology employed for the burst sample 



selection. 

The flrst coarser burst selection was based on the analysis of the GBM telemetry 
packets, which are automatically produced during a trigger and co ntain all trigg er 



20091 ). 



information such as locations, classifications and accumulated rates (iMeegan et al.l 
The so-called GBM "maximum rates" observed over a short period after trigger time 
(<4 s) are produced from the accumulations made for the trigger algorithms and are 
evaluated as statistical significance (signal over noise ratio; SNR) versus the background. 
Typically, the FSW background interval ranges from about —36 s to —4 s with respect to 
the trigger time, thus excluding the most recent few seconds of data and avoiding in most 
cases the contamination by pre-trigger data from the burst. We automatically selected 
bursts showing an increase of more than 80 counts/s over background in at least one BGO 
detector over the full BGO energy range (~250keV to ~40MeV). 

The refined burst selection was based on the analysis of BGO CTIME light curves. As 
previously mentioned, CTIME data have a 64 ms temporal resolution during burst-mode 
and spectral resolution of 8 energy channels. Channel edges are controlled using the specific 
Lookup Tables (LU Ts), which map the 4096 raw channels into the 8 energy channels 



(IMeegan et al. 



20091 ). Exact channel boundaries can vary from detector to detector (BGO 
or 1) and from burst to burst. The BGO CTIME background was computed including 
pre- and post-trigger time intervals, usually from —300 s to +300 s in case of long bursts 
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and from —50 s to +50 s for short bursts, excluding the interval of the burst itself. The 
background-subtracted light curve was then examined for the maximum or "peak" count 
rate on the 64 ms-timescale over each individual CTIME energy channel. 

The resulting total number of GRBs included in this spectral analysis is 52, which 
approximately represents ~20% of all bursts detected during the first year of GBM 
operation. These bursts are listed in Table [TJ The first three columns list the GBM trigger 
number, the GRB namcl, and the burst trigger time (in MET). The numbers of the Nal 
(from to 11) and BGO (0 and 1) detectors used for the temporal and spectral analysis 
are reported in columns 4 and 5. Column 6 gives the angle [9) of the burst with respect to 
the LAT boresight. The LAT FoV covers those events which are located at 6* < 65°. This 
value represents the initial angle from the source calculated at trigger time and can vary 
during the burst in the case of a slew of the spacecraft (the so-called autonomous repoint 
recommendation or ARR). The GBM data type and time interval (with respect to the burst 
trigger time, Tq) adopted for the spectral analysis are listed in the last three columns. 

The full sample of 52 bursts was further subdivided according to the detection 
significance of the event peak in a BGO energy channel. GRBs detected with more than 
3 a significance in the first six BGO energy channels constitute the corresponding Channel 
sample. All bursts in the full sample are detected in Ch.O (~200-500 keV) and Ch.l 
(~500-1000 keV), 28 bursts are detected in Ch.2 (~l-2 MeV), 14 bursts are detected in 
Ch.3 (~2-5 MeV), and 6 bursts are detected in Ch.4 (~5-10 MeV). 

Figure [T] shows an example light curve of one of the brightest GBM bursts in our 



^The naming convention follows the Gamma-ray bursts Coordinates Network (GCN) 
publication policy. Bursts which were not reported in a GCN circular are not given a name 
in column 2. 
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sample, namely the long GRB 081215A, which is detected up to > 5 MeV in BGO. The top 
panel shows the 8-200 keV band, covered by the most illuminated Nal detector (s). The 
other panels show the BGO light curve in different energy ranges, covering five CTIME 
energy channels {Ch.O-Ch.4). This very bright burst was also marginally detected by the 
LAT. Indeed, 9 ~86°, which means that neither directional n or energy info r matio n could 
be obtained with the standard analysis procedures. However, jpelassa et al. ( 2010 ) recently 
presented a new technique to recover the signal from the GRB prompt emission between 
~30 MeV and 100 MeV, which differs from the standard LAT analysis (the so-called "LAT 
Low-Energy" technique, or LLE). Using such non-standard data selection, over 100 events 
above background are detected within a 0.5 s interval in coincid ence with the rnain G BM 



peak. The significance of this excess is greater than 8a (see also 



McEnery et al. 



20081 ). 



3.1. BGO Effective-Area Correction 

In order to correct for the dependence of the BGO effective area on the incidence 
angle, we calculated an additional scaling factor accounting for the angle between the burst 
position and the BGO detectors. This was mainly necessary since the peak count rate 
analysis was perform ed on data without taking the instrument response into account. In 



Bissaldi et al. 



( 120091 ). the off-axis response of the BGO detectors was measured at different 
energies for both flight module detectors at various angles between 0° (i.e., on axis) and 90°. 
A scaling factor could be calculated for each incident direction and then used to correct 
the peak count rates. The correction factors are relatively small out to ~40° and strongly 
increase toward 90°. At higher energies, the correction factor is not as high as at lower 
energies. This mainly reflects the strong absorption of low-energy photons by the BGO 
photomultiplier tubes (PMTs). 



- 10 - 



3.2. LAT Detections 

An interesting result of the BGO peak count-rate analysis emerges by considering only 
those bursts, which are located either inside or at the edge of the LAT FoV. (i.e. 6 <65° or 
65°< 6 < 90° from the LAT boresight, respectively). Figure |5] shows the BGO peak count 
rate measured in Ch.l for 15 and 11 bursts, which respectively fulfill these conditions. The 
GEM trigger IDs and numbers for both subsamples are listed on the top right corner of each 
plot. Circles (green), stars (orange) and squares (red) represent firm, marginal or no LAT 
detections, respectively. The dotted line marks a "detection limit", which was arbitrarily 
placed at 30 (top panel) and 100 cps (bottom panel) in the Ch.l peak count rate. For those 
bursts wit h lower rates in Ch. 1 no detection has been yet reported from the LAT. The very 



promismg 



Pelassa et al. 



(120101 ) LLE technique may recover the signal to confirm the LAT 
marginal detections and even reveal undiscovered emission from some of the Figure [2] BGO 
bursts lying below these arbitrary thresholds. 

The significance of the above analysis lies in its potential to provide a good predictor 
for LAT detections of GEM GRBs. We plan to implement the relevant software into the 
GBM FSW. Thus, starting from the burst location relative to the LAT FoV, the code would 
perform a finer computation of the BGO Ch.l count rate as measured between 500 keV and 
1 MeV. This information would then be sent to the ground, where it would be rapidly and 
automatically analyzed, and subsequently (as the case maybe) provide a prompt alert of 
space- and ground-based GRB observatories. 
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4. Temporal analysis 
4.1. Duration distributions 

We used the CTIME data of the most illuminated Nal detector for each GRB 



to compute the Tgo duration (IKouveliotou et al 



19931) in the BATSE energy range of 



50-300 keV (for comparison reasons). The background model was determined by fitting a 
low-order (< 4) polynomial function over at least a few hundred seconds before and after 
the trigger time. All standard Tgo durations in the 50-300 keV e nergy range were computed 



Mallozzi. Preece. fc Briggs 



with the spectral analysis software RMfilo (version 3.8) (see e. g. 
20051 ). which was also employed for the spectral analysis presented in §5] The resulting 
durations were thus estimated in photon space (i.e., the time it took to collect 90% of the 
burst photons). A different approach (since the BGO data were not yet programmed in 
RMfit at that time) was adopted to compute the burst durations in the 300 keV-10 MeV 
energy range, where we estimated the Tgo's in count space using software developed at 
MPE, thus obtaining a measure of the so-called "BGO-duration" . Figure [3] includes both 
data sets for demonstration purposes. Although these are not directly comparable (photon 
versus count durations) they serve as qualitative trend indicators as we discuss below. 

The initially selected sample of 52 GRBs includes 18 bursts with a duration Tgo< 
3 s (50-300 keV, short GRBs), and 34 bursts with a duration Tgo> 3 s (50-300 keV, 
long GRBs). The three longest bursts in the sample, with durations greater than 100 s, 
are GRB 081009, GRB 090323 and GRB 096018, while the three shortest ones are 
GRB 081226B, GRB 090328B and GRB 090228. The top panel of Figure E] shows the 
distribution of Tgo, calculated in the 50-300 keV and in the 300 keV-lOMeV energy 



^An ad-hoc version of RMFIT for GBM and LAT analysis was developed by the GBM 
Team and is publicly available at http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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ranges, for all 52 bursts. We note that the duration bimodality found in the BATSE data 



( iKouveliotou et al 



19931 ) is evident even in this small sample of GBM data and in both 
count and photon - space durations. We fitted each distribution with a double Gaussian 
function to estimate the medians and deviations. We found that short bursts peak at 
1.2 ±0.3 s (1.04 ±0.16 s) in the 50-300 keV (300 keV-10 MeV) bands, while long bursts 
peak at 33 ±5 s (25 ±8s), respectively. In the bottom panel of Figure [3] we plot the Tgo 
distributions against each other. The dashed line represents a linear fit applied to the data, 
It has a slope of 0.98 ± 0.03, and lies below the bisector (continuous line), indicating that on 
the average, the EGO dura tions are smaller t h an th e Nal's. This result agrees qualitatively 
with the earlier findings of 



Richardson et al. 



( 1l996l ) that GRB durations decline with 



energy. In §4.21 we explore and expand this relationship to the MeV range for the first time. 



4.2. Evolution of Duration with Energy 



We followed the approach described by iRichardson et al.l (Il996[ ) , who presented the 



analysis of 72 intense GRBs from the BATSE 3B catalog. They measured their Tgo in 
four broad energy channels, namely 25-50 keV, 50-100 keV, 100-300 keV, and >300keV. 
The bursts used for their study were uniformly selected by their peak photon flux on the 
64 ms time scale. Thanks to the broader BGO energy coverage, we can measure Tgo over 
the five additional energy channels described in § |3] (the only exception is Ch. 0, which is 
here defined between 300-500 keV, to match the upper edge of the BATSE Tgo values, for 
comparison reasons). The Tgo calculated using Nal only data in the 50-300 keV energy 
band represents our lowest measurement; this energy interval was not further subdivided. 

The energy intervals used for our analysis are listed in the first column of Table 
[21 Columns 2-4 give the number of bursts detected over the different energy intervals, 
where T indicates the total number of GRBs, and L and S indicate long and short bursts. 
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respectively. The average Tgo values calculated using NaT only data in the 50-300 keV for 
each subgroup (again in photon space) are given in columns 5 and 6. It is worth noting 
that the bursts detected up to higher energies are systematically longer. This behavior is 
seen for both long and short GRBs. 

We then fitted the energy versus duration values for each burst in the sample, excluding 
the Nal data points (to ensure a homogeneous set) to a power-law function given by 

T9o = ^o^"^° (1) 

using a minimization technique. We used the central energy value to represent each 
energy channel in the fit. The distribution of the best fit indexes, ago, over all energy 
intervals is shown in Figure |H Blue and red shaded histograms represent the distributions of 
long and short GRBs, respectively. The average value of ago calculated with this technique 
is —0.39 ± 0.04 (0.34), where the first error is the error in the mean, and the error given in 
parenthesis is the standard deviatio n of th e distribution. This value is in good agreement 
with the result of 



iRichardson et al. 



(119961 ) and indicates that the same overall trend of 



declining duration versus energy in the GRB light curves continues up to at least 10 MeV. 

We now proceed to explore the duration-energy relation in subsamples of our 52 GRBs 
selected according to their highest detection in each of the 5 BGO channels described above. 
The sample sizes are defined in Table [2l e.g., 22 long GRBs are detected up to Ch.l and a 
mere 4 up to Ch.5. For each subsample, we first average the values of Tgo for each energy 
channel, and then fit a power law to these data as we did in the total sample. These fits 
were computed separately for long and short bursts and are shown in two panels in Figure |5l 
The values of their ago are given in columns 7 and 8 of Table O The left panel of Figure 
13 clearly exhibits the trend for the longest subsample to extend to higher energies; we can 
see the same trend in the short GRB subsample albeit with lesser statistical significance. 
Whether this trend is simply a selection effect or an intrinsic GRB property remains to 
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be determined with a larger GRB sample, including events with LAT detections possibly 
extending this relation to higher energies. We find evidence for slight curvature of Tqq 
versus Energy in both duration modes indicating a potential cut off. The detailed study of 
these relations using a uniform data set and time-resolved b inning of the NaT du rations is 



201 ll ). We note 



the subject of another publication, currently in preparation (iBissaldi et al.l 
here that the red subsamples include GRBs that have indeed been detected with the LAT 
in GeV energies. 



5. Spectral Analysis 



During the BATSE era, GRB spectra were well represented by a broken power-law 



(the empirical Band function. 



Band et al. 



19931 ) in the BATSE energy band (K06). In 



our spectral analysis, we fitted the time-integrated spectra of the 52 bright bursts in 
Table [T] with two functions: (i) the Band function, and (ii) the Comptonized ( Comp ) 
model. The second represents a low-energy power law with an exponential high-energy 
cutoff, which is equivalent to the Band function without a high-energy power law (i. e., 
(5 — — oo). More details regarding both spectral functions can be found in K06. Since 
we excluded faint or soft GRBs from our sample, the simple power-law model was not 
used. Moreover, no additional extra-components were fitted to both spectral models. For 



a detailed time-integrated spectral analysis of a subsample o 



Guiriec et al. 



thre e short GRBs e x plorin g 



mm . 



Guiriec et al. 



( 1201 ih 



several multi-component emission models see 
have also identified an additional spectral component in the time-integrated analysis of 
GRB 100724B; the appearance of these components in the high-energy part of the GBM 
spectrum could be another predictor for the bursts expected to be seen with the LAT. 

We performed broadband spectroscopy simultaneously on the ~ 8 keV to 1 MeV 
NaT data and the ~ 250 keV to ~38 MeV BGO data. In some cases an effective area 
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correction had to be fit to the BGO data (see e. g. lAbdo et aLl l2009bl ) to match the 
model normahzations given by the Nal data; this correction is usually consistent with the 
uncertainties in the GBM detector responses (<10%). Columns 8 and 9 of Tabled] list the 
time intervals used for each burst. All times are referred to the burst trigger time Tq. 

The time-integrated spectral best fit results are presented in Table [3] for each burst. 
The GRB trigger number is given in column 1, while the best model for the spectral 

umn 

S ..e e.e..e „ea _..o„ ..o.. T.e of . . .... CST^ 

degrees-of-freedom (DOF) is listed in column 9. The Band function was preferred for 
those bursts for which an improvement of > 10 was observed in the CSTAT statistic over 
the Comp function. This assures that the spectra have a well-identifiable high-energy 
power-law component. We find that the Band function is preferred over the Comp model 
in 27 out of 52 cases (52 %). Most short bursts are best fit by a Comp model: Only 5 
out of 18 short GRBs are best f it by the .Band funct ion. Three of those are the brightest 



fit is listed in column 2. The spectral parameters are given in columns 3-7. Co 

i9 



bursts in our sample, for which iGuiriec et al. 



( I2OIOI ) have reported the presence of an 



extra power-law component in the GBM data. We proceed below in the description of the 
spectral parameter distributions, their correlations and their comparisons with empirical 
relationships in the literature. 



5.1. Distributions of Spectral Parameters 



The distribution of the two parameters of the Comp model (index A and -Epcak) are 
shown in Figure El those of the Band function parameters are displayed in Figure [71 In 



^The CSTA T statistics is equivalent to the XSPEC (jArnaudlll996l ) implementation of the 



Cash statistic ((Cash 



19791). 
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both figures the -Epeak parameter distribution is plotted at the bottom. In Figure O the two 
top panels show the low-energy index a and the high-energy index /3. In all plots of the 
following sections, blue, red, and black histograms represent the distributions of 34 long, 18 
short, and the entire sample of 52 bright GRBs, respectively. 

The spectral parameter distributions show that short bursts tend to have larger a (A) 
and smaller /3 values. The four short bursts best fitted with a Band function show higher 
iJpeak than long bursts (Figure [71 bottom panel); in three of them we found a value of 
/3 < —2.6. The Comp model is preferred for the other 13 short bursts. The Comp i?pcak 
distribution of the total sample peaks at ~ 800 keV, while the Band Ep^ak is much lower, 
around 200 keV. Below we compare our results with those of K06. 



5.2. Comparison to the BATSE bright 
burst catalog 

We compare the spectral parameter distributions of our 52 bright GEM bursts with 



the distributions of bursts fro m the BES' 
than the one recently used by iNava et al. 



J sam ple of K06. Our sample is much smaller 



( 120111 ). which comprises the entire GEM GRE 



database. The K06 spectral catalog comprises 350 bright GREs observed with EATSE 
(~20 keV-2 MeV) and is the most comprehensive study of spectral properties of GRE 
prompt emission to date, thus representing a perfect sample for comparing with GEM burst 
properties. Comparison histograms of GEM and EATSE low-energy (a), high-energy {(3) 
and -Epeak spectral parameter distributions are shown in Figure |8l The GEM distributions 
(black empty histograms) are overplotted on the EATSE ones (gray filled histograms). The 



iKaneko et al.l ( 120061 ) designated the model with the more constrained parameters as the 
best-fit (EEST) model. 
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former follow the right y-axis and are rescaled for comparison purposes. 

While the two /3 distributions look similar, differences appear in the a and .Epcak 
distributions. After one year of operations, GBM detected a sample of bright bursts which 
tend to have larger a and higher i?peak values than what was observed in 10 years with the 
BATSE instrument. 

The energy fluence distribution for both samples is calculated in the (BATSE) energy 
range of 25-2000 keV and is also shown in the bottom right panel of Figure El It becomes 
immediately evident that the two samples were selected following different criteria. The 
K06 sample comprises more GRBs with higher fluence than the GBM sample. This 
difference is even more clearly demonstrated in Figure |9l where we plot the GBM burst Tgo 
durations versus their fluences (top right panel) and peak fluxes (calculated over 128-ms, 
bottom right panel) between 8 keV-40 MeV. A comparison of the two panels shows that a 
high-fiuence criterion would have included by far less short GRBs in our sample, as their 
total energies are much smaller than those of the long events, while their peak fluxes span a 
broader range. With 18 short GRBs out of 52 selected bright bursts (~30%), the first year 
of GBM sample contains many more bright short and hard bursts than K06. In fact, only 
17 out of 350 BATSE bright bursts are short, representing ~5% of the sample. Another 
reason for the difference between the two samples is likely the improved trigger algorithms 
implemented in the GBM FSW. BATSE had three time-interval algorithms (based on one 
energy interval), while GBM currently employs 28 trigger algorithms at various time and 
energy channel combinations. These algorithms have vastly improved the capability of the 
instrument to trigger on short and hard GRBs, compared to BATSE. 
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5.3. Correlation among Spectral Parameters — Empirical Relations 

Empirical correlations among spectral parameters have been previously found with 
smaller GRB samples either within individual bursts or for collections of time-resolved 
parameters of multiple events. K06 found no indication of global correlations among the 
time-integrated spectral parameters of their BEST sample (discussed in § I5.2p . while they 
found strong correlations among the time-resolved spectral parameters. They also noted 
that it is best to look for parameter correlations within individual bursts to eliminate 
possible effects due to cosmological redshift that varies from burst to burst. Since no 
time-resolved spectral analysis was performed in this work, we limit our correlation analysis 
to comparisons of the low- and high-energy spectral parameters against i^peak values and 
against each other. 

Figure [10] shows the scatter plot of the Band function high-energy index (3 versus Ep^ak 
for the GEM data alone (top panel), and for the combined data sets of GBM and BATSE 
(bottom panel). The most distinct differences between the two samples is the larger .Epeak 
span of the GBM data, and the larger /3 spread of the BATSE data. The reasons for the 
former have been elaborated in § 15. 2( the latter is the effect again of the K06 selection for 
fiuence and not for hardness (as in the GBM sample) and constitutes, therefore, a more 
representative characteristic of the GRB population as a whole. The GBM subsample of 
the bright, hard events, also includes - not quite unexpectedly - the majority of the LAT 
GRBs, and falls well within the BATSE (3 range. 

Figure [9] displays the -Epeak and Tqq distributions against energy fiuences and peak 
fluxes measured over the entire GBM energy band, i.e., 8keV-40 MeV. -Epeak is plotted 
against fiuences for short and long GBM GRBs (red and blue data points, top left panel) 
and for GBM and BATSE bright GRBs (black and grey data points, bottom left panel). 
Short, low-fiuence bursts show higher -Epeak values, while long, high-fiuence bursts tend to 
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have lower .Epeak values. The distribution of Tgo measured in the standard BATSE energy 
range of 50-300 keV versus energy fiuence in the GBM 8 keV-40 MeV band is shown in the 
top right panel of Figure M This panel clearly exhibits the trend already shown with the 
Swift data that short bursts have lower fluences than longer ones and that high fluences 



unambiguously correspond to longer durations (jPizzichini et al. 



20091 ). The lower right 



panel exhibits that both long and short GRBs have a similar (broad) peak flux range. 

Finally, we plot in Figure [11] the distribution of hardness ratios versus Tgo. The 
hardness ratios are defined by the ratio of counts collected in the BGO energy range over 
those collected in the Nal energy range (1000-40000/8-1000, in keV units). Although the 
current sample is not very large, it allows us to distinguish that shorter bursts in general 
tend to have harder spectra than the long ones. 



Ama ti et al. (2002) 



However, 



Nakar fc Piranl f l2005al ) and 



r eported a re 



ation of E^^^ ve r sus E ^p^k using the BATSE data set. 



Band fc Preecd (120051 ) showed that their results may 



have suffered from strong selection effects and were inconsistent with a larger set of GRB 
data obtained with BATSE. We tested here the Amati relation for those bursts included in 
the BGO-bright burst sample with a reds hift measurement. W e only had 7 cases, namely 



six long GRBs (GRB 080916C, z = 4.35 flGreiner et al.ll2009f ): GRB 090102, z = 1.55; 



GRB 090323, z = 3.57; GRB 090328, z = 0.74; GRB 090424, z = 0.54; and GRB 090618, 
z = 0.54) and one short burst (GRB 090510, z = 0.90). We note that lo n g bur sts nicely 
follow the E'pcak^-E'iso relation, as was also recently pointed out by lAmatil ( 120101 ). and the 
only outlier is the short burst. A recent study by Goldstein et al. (2011) explores the Amati 
relationship in detail using the entire GBM data set. 
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6. Conclusions 

We have studied here a sample of 52 bright and hard GRBs collected during the first 
year of the GBM operation. We have performed temporal and time-integrated spectral 
analysis of all these events and studied the distributions and evolution of the derived 
parameters. The new spectral capabilities afforded with the GBM BGO detectors have 
enabled us to produce a predicting filter using GBM data alone of potential GRB detections 
with the LAT on Fermi. This filter will be implemented in the GBM FSW and alerts will 
be distributed to the wide scientific community to allow timely mult i- wavelength follow up 
observations and, thus, broadband spectral energy distribution studies in GRBs. 

Our temporal evolution analysis has, for the first time, extended the duration-energy 
relationship (found earlier in the BATSE data) to the MeV energy range. Although the 
LAT GeV detections seem to be longer in some G RBs (even up to h undreds of seconds, as 



AbdoetaL 



2011bl ). there seems to be a 



in the case of GRB 090323 and GRB 090328, see 
single power law relation (of index —0.4) between duration and energy in the keV to MeV 
prompt gamma— TSij emission. Whether the GeV emission seen with the LAT in several 
of these GRBs is related to the prompt or the afterglow emission is still an open question, 
which requires more data for definite conclusions. 

Finally, we show that the novel GBM trigger algorithms have improved the collection 
of short and hard GRBs, compared to the BATSE sample. We confirm that their spectral 
parameter distributions are overall similar to those of the K06 sample, and that short 
GRBs are in general harder than longer events. The small subsample of GRBs with known 
distances in our data, is not sufficient to test the various empirical relations in the literature. 
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Fig. 1. — Background-subtracted light curves of the long GRB 081215 observed with the 
GEM detectors. The panels show the sum of the counts in different energy bands as seen 
by {top histogram) the most illuminated NaT detector in the 8-200 keV energy range, and 
{bottom five histograms) , the BOO detectors covering the first five CTIME energy channels. 
The bin width is 128 ms. 
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Fig. 2. — BGO peak count rate measured in Ch.l (500 keV - IMeV) for those bursts inside 
{top panel) and at the edge {bottom panel) of the LAT FoV. Circles (green), stars (orange) 
and squares (red) represent firm, marginal or no LAT detections, respectively. The dotted 
line marks the arbitrary "detection limit" placed at 30 (top panel) and 100 (bottom panel) 
count s/s. 
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Fig. 3. — Top Panel: Burst duration distributions in the 50-300 keV {gray filled histogram) 
and in the 300 keV-lOMeV {black empty histogram) energy ranges for 52 bursts of our 
bright-burst sample. Both distributions show a bimodal shape. Bottom Panel: Scatter plot 
of the burst duration distributions. The continuous line represents perfect linearity, while 
the dotted line represents the fit applied to the data. 
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Fig. 4. — Distribution of the best fit indexes ago of power-law fits to the burst durations 
over different energy channels. Blue and red histograms represent the distributions of long 
and short GRBs, respectively. The black histogram represents the entire sample. 
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Fig. 5. — Evolution of the Tgo duration with respect to energy for four subsamples of 12 
long GRBs (left panel) and 7 short GRBs (right panel). Energy bands are over the 5 BGO 
CTIME channels. Different curves correspond to different burst sets for which such duration 
measurements were possible (see also text and Table [2]). 
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Fig. 6. — Index A {top panel) and Ep^ak {bottom panel) distribution of the time-integrated 
spectra which are best fitted with the Comp modeL Blue and red histograms represent the 
distributions of long and short GRBs, respectively. 



- 31 - 



BAND 




Fig. 7. — Low-energy index a, high-energy index /? (top panels) and -Epeak (bottom panel) 
distribution of the time-integrated spectra best fitted with the Band function. Blue and red 
histograms represent the distributions of long and short GRBs, respectively. 
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Fig. 8. — Low-energy index a (top left panel), high-energy index j3 (bottom left panel), -Epcak 
(top right panel) and energy fluence (bottom right panel) distributions of time-integrated 
spectra from 350 bright BATSE bursts (grey filled histogram) from K06 (BEST sample) 
and 52 bright GBM bursts (black histogram). The GBM energy fluence is computed in 
the standard BATSE energy range, namely between 25 and 2000 keV. The GBM parameter 
distributions follow the right y-axis. 
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Fig. 9. — Left panels: Distribution of .Epcak versus energy fluence calculated over the en- 
tire GBM energy range (top left panel, 8keV-40MeV) and over the BATSE energy range 
(bottom left panel, 25-2000 keV). Blue and red data points represent long and short GRBs, 
respectively; Grey data points represent 350 GRBs from the K06 sample. Right panels: 
Distribution of Tgg (50-300 keV) versus energy fluence (8 keV-40 MeV, top right panel) and 
versus 128 ms peak photon fluxes (8keV-40MeV, bottom right panel). 
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Fig. 10. — Scatter plot of the GBM spectral parameter pair E^cak^P for 22 long GRBs and 
5 short GRBs (blue and red data points in the top panel). Gray data points in the bottom 
panel represent 350 GRBs from the K06 sample. Bursts fully or marginally detected by the 
LAT are marked with green stars. 
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Fig. 11. — Distribution of hardness ratios versus Tqq (50-300 keV). The hardness ratios 
are defined as the ratio of counts collected between 1-40 MeV over those collected between 
8-1000 keV. Bursts fully or marginally detected by the LAT are marked with green stars. 
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Tablc 1. Basic properties of 52 bright GRBs 
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Note. — ^ The time range values are given in s relative to the trigger time Tq. They represent 
the interval used for the time-integrated spectral analysis 
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Table 2. Duration versus Energy parameters 



Energy 
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300 keV- 


- 2 MeV 


9 6 


3 


43.1 ±0.7 


1.08 ±0.18 


-0.41 ±0.04 - 


0.42 ±0.17 


300 keV- 


- 5 MeV 


4 2 


2 


52.9 ±0.5 


0.92 ±0.15 


-0.38 ±0.04 - 


■0.42 ±0.17 


300 keV- 


10 MeV 


6 4 


2 


76.7 ±0.6 


2.24 ±0.16 


-0.32 ±0.02 - 


0.40 ±0.12 



Table 3. Summary of time-integrated spectral fit results OF 52 bright GRBs 



GBM 


Best 


A 


-^pcak 


index 






AcflF 


CSTAT/ 


TVig. # 


Model 


(ph. cm~^ keV~^) 


(keV) 


Corr. 


DOF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


080723.557 


Band 


0.035 ± 0.0008 


219 (tl) 




-0.858 (l°-°i) 


-2.63 (l°-i^) 


0.83 


384/235 


080723.985 


Comp 


0.0111 ± 0.0003 


445 {tli) 


-0.949 (10-°i) 




496/359 


080725.541 


Comp 


0.0135 ± 0.0009 


1670 (+Ho) 


-0.89 (+0-}^) 








362/359 


080802.386 


Comp 


0.0166 ± 0.0025 


600 {t\To) 


-0.65 (t°i^) 








402/358 


080807.993 


Comp 


0.0043 ± 0.0002 


790 


-1-01 n'oi) 








577/481 


080816.989 


Comp 


0.0039 ± 0.0002 


1550 


-0.51 (;S:n) 








390/360 


080817.161 


Band 


0.0147 ± 0.0004 


425 




-0.99 ± 0.02 


-2.3l'(i[]-i^) 

V U--L4:/ 




544/361 


080825.593 


Band 


0.0617 ± 0.0033 


176 it 7 




-0.64 ± 0.04 


-2.52 (+°-}°) 

^ V— 0.14/ 


0.80 


418/357 


080905.499 


Comp 


0.0090 ± 0.0016 


500 (+]^?^ 

<juu V — 110/ 


-0.20'(«j9) 






471/361 


080906.212 


Band 


0.1152 ± 0.0110 


147 (+^) 




-0.37 ± 0.06 


-2.39 (+°-l°) 

V— 0.12/ 




522/479 


080916.009 


Band 


0.0166 ± 0.0002 


540 f+32\ 




-1.06 ± 0.02 


-2 24 (+°-?^) 
^•^^ V-O.IO/ 




533/358 


080925 775 


Band 


0272 ± 0016 


138 ± 8 




-0 96 lb 04 


-2 27 (+°-°^) 

^,-0.10^ 




430 /359 


081006.604 


Comp 


0.0038 ± 0.0004 


840 (+^^0) 


-0.43 (+0-3 ) 






423/362 


081009.690 


Band 


0.0296 ± 0.0143 


63 (1}^) 




-0.53' (l°:L) 




0.70 


502/358 


081012.045 


Comp 


0.0123 ± 0.0011 


750 itl^) 


-0.44' (i°:i) 








368/355 


081024.891 


Comp 


0.0096 ± 0.0014 


1300 


-0.46 (1°:,\) 








326/359 


081101.532 


Comp 


0.0251 ± 0.0008 


290 ± 30 


-0.686 ± 0.04 








360/360 


081110.601 


Comp 


0.0105 ± 0.0004 


470 


-1.064 (t^:"!!) 


-0.47' (l^i^) 






410/359 


081121.858 


Band 


0.0266 ± 0.0049 


158 (ii?) 


-1-94'CT 




456/356 


081122.520 


Comp 


0.0093 ± 0.0006 


207 {tlD 


-0.90 ± 0.07 








435/360 


081125.496 


Band 


0.0932 ± 0.0051 


164 ± 7 




-0.48 ± 0.05 


-2.86' (l°i^) 


0.76 


408/358 


081126.899 


Comp 


0.0040 ± 0.0004 


330 (11°) 


-1.00 ± 0.07 








488/360 



Table 3 — Continued 



GBM 


Best 


A 


keV j 


-^peak 


index 




a 






Acff 


CSTAT/ 


TV,^ ^ -LL 

ii-ig. # 


Modei 


(ph. cm" 


(keV j 






L^orr. 


JJUr 


(V) 




f3) 


y*) 




(6) 


(7) 




f9) 


081129.161 


Band 


0.0114 ± 


0.0007 


250 ± 30 




-0.96 


V-o.osJ 


-2.20 


/'+0.15\ 
1-0.22^ 




454/356 


081207.680 


Band 


0.0095 ± 


0.0002 


430 {tf^ 




-0.67 


± 0.03 


-2.11 


/+0.08\ 
1-0.09^ 




501/356 


081209.981 


Band 


0.0347 ± 


0.0021 


1100 (1^0°) 




-0.68 


(+0.1A\ 

\— U. ii/ 


-2.23 


(+Q.2&\ 
V— U.o / 




349/357 


081215.784 


Band 


0.1074 ± 


0.0014 


458 ± 13 




-0.71 


/+0.02\ 
I n no I 


-2.36 


(+0.0A\ 
\— U.UO/ 




519/356 


081216.531 


Band 


0.0200 ± 


0.0009 


1270 [t^Z) 




-0.81 


/'+0.07\ 
\ — U.Uo/ 


-2.54 


/'-l-0.40\ 
V-O.80J 




453/357 


081224.887 


Comp 


0.0378 ± 


0.0007 


380 ± 11 


-0.73 ± 0.02 










0.87 


442/358 


081226.509 


Comp 


0.0311 ± 


0.0044 


350 +^°° 


-0-41 (l|]i^) 








;^_0.50/ 




414/360 


081231.140 


Band 


0.0149 ± 


0.0006 


234 (1^6) 




-1.06 


± 0.04 


-2.26 




472/358 


090102.122 


Comp 


0.0180 ± 


0.0003 


412 (li^) 


-0.86 ± 0.02 










0.87 


442/357 


090131.090 


Band 


0.0321 ± 


0.0030 


61 ± 4 




-1.21 


/+0.08N 
0.06/ 


-2.42 


/+0.09N 
^—0.10/ 




348/236 


090217.206 


Comp 


0.0125 ± 


0.0002 


633 11) 


-0.91 ± 0.02 












659/479 


090219.074 


Comp 


0.0289 ± 


0.0115 


^ , /_Li nnx 

214 i^n 


„ /_i_n s\ 

-0-2 [t'oD 












353/361 


090227.310 


Comp 


0.0055 ± 


0.0002 


1000 (1?00) 


-0.86 ± 0.06 












383/360 


090227.772 


Band 


0.0762 ± 


0.0016 


2100 ± 100 




-0.51 


/+0.03N 
1-0.02^ 


-3.33 


('+0.27N 
I-O.40J 




548/479 


UyU2zo.2U4 


Comp 


O.Oroo ± 


0.0016 


o4(J ± oU 


-U.biJ ± U.Uo 








r Ac\ / /ion 

O4U/450 


090305.052 


Comp 


0.0126 ± 


0.0006 


890 (11^°) 


-0.58 (l°:J^) 












570/480 


090308.734 


Comp 


0.0193 ± 


0.0080 


664 (1^°) 


-0.53 (IH?) 








/'+o.09^ 
I-0.12; 




688/600 


090323.002 


Band 


0.0178 ± 


0.0003 


530 (li) 




-0.81 


± 0.02 


-2.42 




568/237 


090328.401 


Band 


0.0173 ± 


0.0003 


660 ± 40 




-0.93 


± 0.02 


-2.44 


/'+0.14\ 

1-0.19; 




534/360 


090328.713 


Comp 


0.0319 ± 


0.0016 


2000 (li°) 


-0.96 ± 0.07 












376/356 


090330.279 


Band 


0.0068 ± 


0.0005 


198 ± 19 




-0.92 


/'+0.07^ 
l.-o.oe/' 


-2.28 


/'+0.15N 
U0.21; 




747/477 


090424.592 


Band 


0.1419 ± 


0.0026 


146.2 (11^) 




-0.86 


± 0.02 


-2.76 


± 0.08 


0.78 


843/478 


090429.753 


Comp 


0.0150 ± 


0.0008 


1400 (112°) 


-i-06'(lS:g^) 












372/360 



Table 3 — Continued 



GBM 
Trig. # 

(1) 


Best 
Model 
(2) 


A 

(ph. cm~^ keV~^) 
(3) 


-^peak 

(keV) 


index 
(5) 


a 
(6) 


(7) 


Aeff 

Corr. 
(8) 


CSTAT/ 
DOF 
(9) 


090510.016 


Band 


0.0427 ± 0.0009 


4400 (tTo) 




-0.79 ± 0.03 


-2.74 




486/479 


090528.516 


Band 


0.0170 ± 0.0070 


170 {tlD 




-1.10 ± 0.04 


-2.21 (in^) 




553/358 


090531.775 
090618.353 


Comp 
Band 


0.0094 ± 0.0006 
0.0717 ± 0.0011 


1750 

150.2 ± 2.7 


-o.m'itrol) 


-1.12 ± 0.01 


-2.50 ± 0.03 


0.78 


632/481 
532/233 


090620.400 


Band 


0.0920 ± 0.0045 


150 ± 4 




-0.174 ± 0.05 


-2.77 it^il) 


0.70 


580/479 


090623.107 


Band 


0.0084 ± 0.0004 


307 {tlD 




-0.63 (^HD 


-2.29 


0.72 


607/478 


090626.189 


Band 


0.0440 ± 0.0014 


160 (+^) 




-1.04 ± 0.03 


-2.30 


0.74 


353/240 



Note. — ^ The time range values are given in s relative to the trigger time Tq. They represent the interval used for 
the time-integrated spectral analysis 



